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Cylindrical spacetimes with rotation are studied using the Newmann-Penrose formulas. By study- 
ing null geodesic deviations the physical meaning of each component of the Riemann tensor is given. 
These spacetimes are further extended to include rotating dynamic shells, and the general expression 
C^ , of the surface energy-momentum tensor of the shells is given in terms of the discontinuation of the 

f^ ' first derivatives of the metric coefficients. As an application of the developed formulas, a stationary 

f^ I shell that generates the Lewis solutions, which represent the most general vacuum cylindrical solu- 

04 ■ tions of the Einstein field equations with rotation, is studied by assuming that the spacetime inside 

,— ( ' the shell is flat. It is shown that the shell can satisfy all the energy conditions by properly choosing 

^^ , the parameters appearing in the model, provided that < cr < 1, where a is related to the mass per 

' • unit length of the shell. 

00 ■ PACS numbers: 04.20Cv, OA.SO.+x, 97.60.Sm, 97.60.Lf. 

cni : 

^ ; I. INTRODUCTION 

00 ; 

C^^ . Gravitational collapse of a realistic body has been one of the most thorny and important problems in Einstein's 

^^ ' theory of General Relativity. Due to the complexity of the Einstein field equations, the problem even in simple cases, 

such as, spacetimes with spherical symmetry, is still not well understood M, and new phenomena keep emerging 

0. In 1991, Shapiro and Teukolsky |g] studied numerically the problem of a dust spheroid, and found that only 

^^ the spheroid is compact enough, a black hole can be formed. Otherwise, the collapse most likely ends with a naked 

■ — «. I singularity. Since then, the gravitational collapse with non-spherical symmetry has attracted much attention. In 

pl^' particular, by studying the collapse of a cylindrical shell that consists of counter-rotating particles, Apostolatos and 

I Thorne (AT) found that the rotation always halts the collapse [Q. As a result, no naked singularities can be formed 

wp.' on the symmetry axis. However, in the AT work only the case where the total angular momentum of the collapsing 

. . , shell is zero was considered. In more realistic case, the spacetime has neither cylindrical symmetry nor zero angular 

momentum. As a generalization of the AT work, in this paper we shall consider the case where the total angular 

momentum is not zero, while still keep the requirement that the spacetime be cylindrical. 

Another motivation for us to consider rotating shell comes from recent study of the physical interpretation of the 
d ' Levi-Civita vacuum solutions H, and the Levi-Civita solutions with cosmological constant (LCC) M. These solutions 
have been known for a long time |Q], but their physical properties were studied extensively only very recently [^, 
and there are still several open problems to be solved. By looking for some physical sources to the LC solutions, we 
extended one of the two parameter that is related to, but not equal to, the mass per unit length, from the range, [0, j], 
to the range, [0, 1] M. However, beyond this range, the physical meaning of the solutions is still not clear. Despite of 
the simplicity of the solutions, it was found that they have very rich physical meaning. In particular, they can give 
rise to black hole structures 0. Since the LC and LCC solutions are all static, it is very interesting to generalize 
these studies to the rotating case. 

As it can be seen from the discussions given below, when rotation is included, the problem is considerably compli- 
cated. This partially explains why spacetimes with rotation are hardly studied (analytically). Thus, to start with, 
in the next section (Sec. II) we shall first study the main properties of cylindrical spacetimes with rotations, using 
the Newmann-Penrose (NP) formulas M. One of the main reasons to use the NP formulas is that they give directly 
physical interpretation for each component of the Riemann tensor. This is particularly useful when spacetimes contain 
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gravitational waves. In Sec. Ill, the spacetimes are extended to include the case where rotating matter shells appear. 
To deal with the problem, one usually uses Israel's method [Q. However, we find that for the present problem Israel's 
method becomes very complicated and is very difficult to be implemented. Instead, we shall follow Darmois ||ll| and 
Lichnerwic z II4 (see also Papapetrou and Hamoui [Q and Taub [Q). Although the two methods are essentially 
equivalent [|15[, the latter is simpler, specially in dealing with complicate boundaries like the present one. The dis- 
advantage of the latter is that it requires only one set of coordinates across the boundaries, while Israel's method 
does not. In this section, the physical interpretation of the surface EMT is also studied by solving the corresponding 
eigenvalue problem. As an application of the developed formulas, in Sec. IV we consider a stationary shell that 
generates the Lewis solutions, which represent the most general cylindrical vacuum solutions of the Einstein field 
equations with rotation |16|| , while in Sec. V our main conclusions are given. 



II. CYLINDRICAL SPACETIMES WITH ROTATION 



To begin with, let us consider the cylindrical spacetimes with rotation described by the metric |T^, 

ds^ = e^^X-'A) (rfi2 _ ^^2) _ ^2g-2^ (^^^ _^ ^^)2 _ ^2^^^2^ 



(1) 



where ip, Xy ^ ^-i^d uj are functions t and r, and {a;''} = {t,r, z,ip}, (fi — 0,1,2,3), are the usual cylindrical 
coordinates. In general the spacetimes have two Killing vectors, one is associated with the invariant translations 
along the symmetry axis, ^(2) = dz, and the other is associated with the invariant rotations about the axis, ^(^) = dip. 
Clearly, for the metric given above, the two Killing vectors are orthogonal. When lu ^ 0, the metric represents 
spacetimes without rotation, in which the polarization of gravitational waves has only one degree of freedom and the 
direction of polarization is fixed |lq,19|. For the spacetimes to be cylindrical, several criteria have to be satisfied [ pO[ . 
When the symmetry axis is regular, those conditions are easily imposed. However, when it is singular, it is still not 

clear which conditions should be imposed |^l|]. 

Co rres ponding to the metric (|^) the Christoffel symbols and Einstein tensor are given, respectively, by Eqs.( |A.2] ) 
and (A^). To further study the spacetimes, let us consider the Ricci and Weyl scalars, as these quantities have their 
explicit physical interpretations H , and they are particularly useful in the study of gravitational interaction among 
mater fields and gravitational waves 1 22 2^ . Choosing the null tetrad |j 



,..— iu^.o, -.}.„- ^-^ 

1 ( , e"^ 



{1, -1, 0, -iu}.. 



1 



ai' 



^^ - -^ <! e"'^ -i^— 
V2r' ' ' W 



(2) 



where a bar deno tes t he complex conjugate, and A is an arbitrary function of t and r, we find that the spin coefficients 
are given by Eq.(BJ). Since k, — v = 0, the two null vectors l^ and n'* define two null geodesic congruences |^^. The 
one defined by Ijj, is outgoing, while the one defined by n^ is ingoing, similar to the no- rotating case [M. In fact, it 
is easy to show that 
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where e and 7 are given by Eq.(B.l). Choosing the function A as ^ = e-^ ''', from Eq.(B.l) we find that e = 0, and 



then Eq.(|3|) shows that the null geodesic congruence defined by Z^ is affinely parameterized. Hence, the corresponding 
expansion, rotation and shear of the outgoing null geodesic congruence are given, respectively, by 
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^It should be noted that the choice of the null tetrad used here is different from the one used in |24| for the case to 
can be obtained one from the other by exchanging the roles of the two null vectors, /^ and n^. 
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and W^t = dW/dt, W^r = dW/dr, etc. Similarly, if we choose the arbitrary function A as A = e^ ^, Eq.(BJ) yields 



7 = 0, then Eq.fl) shows that the null geodesic congruence defined by n^ now is affinely parameterized, and the 
corresponding expansion, rotation and shear of the ingoing null geodesic congruence are given, respectively, by 



Re{^i) 



a'^d'-X) \Y 



W_r 



2V2 



:,n^^;.)n'''''' - el 



1/2 



W 



.[^^.P^^"^ - [Im{^,)Y - 0, 



2%/2 



aiV-,* - V'.r) - 



w, - W, 



w 



(6) 



Once the spin coef ficie nts ar e giv en, we can calculate the corresponding Ricci and Weyl scalars, which are given, 
respectively, by Eqs.( |B.2| ) and (B.3). From those expressions we can see that all the Weyl scalars are non-zero, and 
each of them has the following physical interpretation |E3,pa: The ^0 and ^1 terms represent, respectively, the 
transverse and longitudinal gravitational wave components along the null geodesic congruence defined by n^, and the 
\l/4 and ^3 terms represent, respectively, the transverse and longitudinal gravitational wave components along the 
null geodesic congruence defined by l^ , while the ^2 term represents the "Coulomb" component. 

The physical meaning of the Weyl and Ricci scalars can be further studied from geodesic deviations. Because of 
the symmetry, it is sufficient only to consider the null geodesies defined by l^ , which are affinely parametrized when 
A — e^^^, as sh own above. Let 77^ be the deviation vector between two neighbor geodesies, and rj'^lfi — 0. Then, 



using Eqs.(B.2) -(B.7), we find that the geodesic deviation can be written in the form 









{$ooe; 







*oe^'' +*(*!+ $01) er 



'03 -ri(*i + $oi)e5'3}?7,., 
where ^q, *i, $00, and <i>oi are given by Eqs. (|B.3[) and (|B.4[ ) with A = e^~^ , and 
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Let So and Sp be infinitesimal 2-elements spanned by ej 



Eqs.(R) has the following physical interpretation p2, 
and 63 and orthogonal to a null geodesic C defined by ^^, passing 5*0 and Sp at the points O and P, respectively. 
Let S be an infinitesimal circle with center O, lying in So as illustrated by Fig.l. Suppose that a light beam meets 
So in the circle 5*, then each of them has the following effect on the image of the circle S on Sp. The first term 
<I>oo in Eq.(0) will always make the circle contracted, as for any physically realistic matter field we have $00 > [cf. 
Fig. 2(a)]. The second term \l/o will make the circle elliptic with the main major axis along ei^, as shown by Fig. 2(b). 
To see the physical interpretation of the last two terms, let us consider a tube along the null geodesic C. Consider a 
sphere consisting photons, which will cut So in the circle S with the point O as its center, as shown in Fig. 3. Then, 
the last term in Eq.(|j) will make the image of the sphere at the point P as a spheroid with the main major axis 
along a line at 45" with respect to e^' in the plane spanned by e'^ and Cg, while the rays are left undeflected in the 
63 -direction. This can be seen clearly by performing a rotation in the plane spanned by e^' and 63 , 



e^ = cos ae I — sm ae 3 , e^ — sm ae ^ + cos ae 
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which leads to 



ei'g = siii(2a)(e'5'e'i - e%e%) + cos{2a)ie"^e% + e'^e'"). (11) 

Thus, choosing a — 7r/4, we have 

ei3 =6161 -6363, (a = 7r/4). (12) 

Combining Eq.([7|) and the above equation, we can see that the last term will make a circle in the e'ie'3-plane into 
an ellipse with its main major axis along the e'l-axis, which is at 45*^ with respect to the ei-axis. It should be noted 
that in the case of timelike geodesies the \l/i term deflects the sphere into an ellipsoid |2^. Moreover, the third term 
in Eq. (|7|) is absent in the timelike geodesic case. The effect of this term will make a clock "flying" with the photons 
slow down, in addition to the effect of deflecting the photons in the e3-direction. It is interesting to note that there 
is a fundamental difference between the time delay caused by this term and the one caused by a Lorentz boost. The 
latter, in particular, has no contribution to geodesic deviations, timelikc or null. From the above analysis we can see 
that for a pure Petrov type N gravitational wave propagating along the null geodesic congruence, in which only the 
component ^p is different from zero, the gravitational wave has only one polarization state, similar to the case without 
rotation [|l8| , P2|J24| • The difference between these two cases is that in the case without rotation, the polarization angle 
remains the same even in different points along the wave path, while when w 7^ 0, in general this is no longer true. In 
fact, it is easy to show that 

e'2-J'' = 0, e^J'^ = ^We^-^^ (e(,' + e^) . (13) 

Thus, although 63 is parallel-transported along the null geodesic congruence, 63 in general is not, and is rotating with 
respect to a parallel-transported basis. Since the polarization angle of the ^g wave remains the same with respect to 
63 , the polarization direction is also rotating with respect to the parallel-transported basis. 

III. ROTATING CYLINDRICAL SHELLS 

In the last section, the main properties of the spacetimes with rotation have been studied. Since the Einstein field 
equations are all involved with derivatives of the metric coefficients up to the second order, and the Bianchi identities, 
which are usually considered as representing the interaction among gravitational fields and matter fields, up to the 
third order, it is generally assumed that the metric coefficients are at least C'^ [p^, that is, the derivatives of the 
metric coefficients exist and continuous at least up to their third order. However, for a long time it has been realized 
that this condition is too strict and rules out many physically interesting cases, such as, shells, star, and so on. 

In this section, we shall generalize the formulas given in the last section to the case where the metric coefficients 
are C^ only in certain regions, while across the boundaries that separate these regions they are C°, that is, the metric 
coefficients are only continuous across the boundaries. These boundaries can be classified into two different kinds, one 
is boundary surfaces, like a star, the other is surface layers, like a matter shell |10|| . In the former case, the extrinsic 
curvatures of the boundaries in their two faces are equal, while in the latter case they are not, and as a result matter 
shells in general appear on these boundaries. In this paper, we shall treat the two cases together and consider the 
former is a particular case of the latter. 

In a given spacetime, there may exist many disconnected boundaries. However, in the following we shall consider the 
case where there exists only one boundary, as its generalization to the cases of many boundaries is straightforwards. 
Assume that the whole spacetime is divided into two regions V by a hypersurface S, where 

y+ = {x^ : > 0}, y^ = {x^ : (/)< 0}, S == {x^ : == 0}, (14) 

with 

</. = r-i?(i), (15) 

where R{t) is an arbitrary function that describes the history of the boundary. Then, for any of the metric coefficients, 
which is C" across the boundary and C^ in the regions V^ can be written in the form pq | 

/(i,r) = /+(i,r)iJ(0) + /-(i,r)[l - i/(0)], (16) 

where / = {i/i, X: W^j ^}; H{(/)) is the Heavside function defined as fj 



^It should be noted that the exact value of H{(p) at the point (f> — is not uniquely defined and can be given any value. 
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and /+ (/ ) is the function defined in the region V^ {V ), with the C° condition 

Hni f^{t,r)= hni f^{t,r). 

Using the distribution theory, it can be shown that 
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where an over-dot denotes the ordinary differentiation with respect to t, (5((/)) the Dirac delta function, and 
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It should be noted that in deriving Eg. ([l9|) we have used the relation 

[f,t]-=-R{t)[fA-. 
Substituting Eq.(|l9|) into Eq.( [A.3| ), we find that the Einstein tensor in general can be written in the form 
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where G+ {Gni,) is the Einstein tensor calculated in the region y+ {V ), and ^^v is the Einstein tensor calculated 
on the hypersurface r — R{t). The non- vanishing components of 7^^ in the present case are given by 
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where the quantities with the subscript "0" denote the ones calculated on the hypersurface r — R{t), for example, 
Wo = w(t, R{t)), and so on. Writing the energy-momentum tensor (EMT) in a form of Eq.(p2|), we find that the 
Einstein field equations G^i^ — kT^^^, where k is the Einstein constant, can be written as 



G+. - A;T+ , (r > i?(i)), 
liLv = kT^^y, (r = i?(i)), 



(24) 
(25) 
(26) 



where r^^ can be interpreted as representing the surface energy-momentum tensor, which in the present case takes 
the form 



where 
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p, ^ (1 - 7i'2)e2(V>o-xo) |[^^^]- _ 2[V.,.]- + K 
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^^ = e'^°J^, (^^ = l^oe-^''{a;o^*+5^}, (29) 



uxu^ = -zxz^ = -Lpx(p^ = 1, 
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with the properties 



In order to have th e phys ical interpretation for each term appearing in Eq.(p7D, we need to cast the surface EMT 
in its canonical form [ P7y2§[| , that is, we need to solve the eigenvalue problem, 

T.^r = A^^ (31) 

This system of equations will possess nontrivial solutions only when the determinant det\T^ — A(5^| = 0, which in the 
present case can be written as 

A(p, - A) [A2 - (,y - p^)\ + q^- r,pA = 0. (32) 

Clearly, the above equation has four roots, A = 0, Pz, A±, where 

1 [, 



'--2 



(r;-p^)±i?i/2j, D={T^+p^f-Aq\ (33) 

It can be shown that the eigenvalue A = corresponds to the eigenvector ^^ = 77'^, where 77^ is the normal vector 
to the hypersurface r = R{t), and given by Eq. (p9|) . The eigenvalue \ = Pz corresponds to the eigenvector ^^ = z'^, 
which represents the pressure of the shell in the z-direction. 

On the other hand, substituting Eq.(|33|) into Eq.(pl[), we find that the corresponding eigenvectors are given, 
respectively, by 

ei = {^±+Pv)u^ + q^''- (34) 

To further study the physical meaning of A±, it is found convenient to distinguish the three cases: (a) D > 0; (b) 
D = 0; and (c) D <0. 

Case (a): In this case, the two roots A± and the two corresponding eigenvectors ^^ are all real and satisfy the 
relations, 

(A+ +p^){\^ +p^) = q^, 

i5i/2(A±+p^) " ' 

e+£.-9^.. = 0. (35) 

From these expressions we can see that when A+ +p^p > 0, the eigenvector ^^ is timclike, while 1^^ is spacelike. Setting 
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we find that Ei^,, {a — 0, 1, 2, 3) form an orthogonal base, i.e., E}^.E(i,^x = ^067 where 
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Then, in terms of these unit vectors, the surface EMT given by Eq.(27) takes the form 
where 

- ^ ^^- +''^^- [b'I-'p^ + {p^{^ + P^) - 2g2] } , (A+ + p^ > 0). 
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Hence, in terms of its tetrad components, t^^ can be cast in the form. 
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This corresponds to the Type I fluid defined in ^^. Thus, in this case the surface EMT represents a fluid with its 
surface energy density given by S, measured by observers whose four- velocity are given by i?/^-,, and the principal 
pressures in the directions E^^. and E^^. , given respectively by p^ and p(3) . 

When A+ +pi^ < 0, the eigenvector ^',1 is spacelike, while ^^ is timelike. Now if we define the two unit vectors i?(o) 



and i<^(3) as 
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we find that the surface EMT also takes the form of Eq. (Eoh , but now with 
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Case (b): In this case we have 
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Then, the two roots A± degenerate into one. It can be shown that this multiple root corresponds to two independent 
eigenvectors, 
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which are all null. From these two null vectors we can construct two unit vectors, one is timelike and the other is 
spacelike. But, these two unit vectors are exactly u^ and <p^. Then, in the base E^^. = {u^, r]^ , z^, 'P^}^ the surface 
EMT wiU take the form 
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with 
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This corresponds to the Type II fluid defined in [E^ . 

Case (c): In this case the two roots X± are complex, and satisfy the relations A_ = A-|_. The two corresponding 
eigenvectors, given by Eq.(p3), now are also complex. This means that in the present case the surface EMT cannot 
be diagonalized (by real similarity transformations), and is already in its canonical form. Thus, in the base consisting 
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to ' (a) 
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from which we can see that now 77 denotes the surface energy density of the shell, measured by observers whose 
four- velocity are given by u^, pz and p^p the principal pressures in the directions, z^ and (p'^ , respectively, and q the 
heat flow in the i^s'^-direction. 

It should be noted that all the above physical interpretations are valid, provided that the surface EMT satisfies 
some energy conditions | |27[ |. 

IV. STATIONARILY ROTATING SHELLS 

As an application of the formulas developed in the last section, in this section we shall consider stationarily rotating 
cylindrical shells. It should be noted that such shells were studied previously by various authors ||29| . However, an 
important difference in the present case is that now we allow the spacetime inside the shell be rotating. 
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where fl is a constant that represents the angular velocity of the uniformly rotating coordinate system pG] . Note that 
in the above metric the axis r = is well-defined and free of any kind of singularities and sources. When r > r2~^, 
the metric coefficient g^Q becomes negative and the Killing vector ^/^^ = 5q becomes space-like. In the following we 

shall assume that the above metric is valid only for r < il~^. 

If we further assume that outside of the shell, the spacetime is vacuum, then it should be described by the most 
general vacuum Lewis solutions of the Einstein field equations, given by p| 



dsi 



fdr + 2kdtdif - h{df^ + dz^) - ld(p^ 



(49) 



where 



/ ^ afi-« _ ^fi+", k^Af, /i = f("'-i)/2, 



/ 



A^f; A 



pl+n 



anf 



(50) 



with a, 6, c and n being constants that can be real or complex. When they are complex, certain relations have to 



For the present purpose, we shall 



:ney are 
be satisfied among them, in order to have the metric coefficients be real |l^j3l| ]. 
consider only the case where they are all real. 

It should be noted that the above expressions are valid only for n 7^ 0. However, the solutions for n = can be 
obtained from them by first letting c = nc and then taking the limit n ^ 0. Since this process is straightforwards, in 
the following we shall consider only the solutions given by Eqs.(E3) and (pw, and consider the ones for n = as their 
particular case. The physical meaning of the parameters a, b, c and n were first studied by Lewis ||16| , and more 



recently by da Silva et al. |3^ . When all the parameters are real, it can be shown that the spacetime is asymptotically 
flat as f — > +00, and singular on the hypersurface f = 0. This can be seen, for example, from the Kretschmann scalar, 

R^..p^R^ = -^;^^ . (51) 

Thus, the Kretschmann scalar diverges as r ^ 0, except for the cases n = ±1. In the last two cases, it can be shown 
that the spacetimes are flat and belong to Petrov type O. When n = 0, 3, the spacetimes are Petrov type D, while 
all the rest are Petrov type /. 

In order to apply the formulas developed in the last section, we need first to write the Lewis solutions in the form 
of Eq.ffl). To this end, we make the following coordinate transformations 
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for n ^ —1, where a, /?, 7, S and d are arbitrary constants, and Rq = [a(l + n)^/4]'*/(^+"^ . Then, the Lewis solutions 
takes the form, 
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dsl = a^A^''{r){dt^ - dr^) - P^A^^<-^^-^\r)dz^ - ^A^'-^-^^\r){dip + 5dtf , {a ^ 1/2), (53) 



where 



^(r) = (r + d)l/(2-l)^ a = aRl'', a=i-^. (54) 



From the above expressions we can see that the range of f, f e [0, 00), is mapped into the range r G [— d, 00). 
Hence, when r -^ +00, the spacetime is asymptotically flat, and when r -^ —d it is singular, except for the case n = 1 
where the spacetime is flat. As we shall see blow, we will use the above metric as describing the spacetime outside 
of a stationary shell. Thus, to prevent spacetime singularity from happening outside the shell, we shall require that 
r + d> 0. 

When a = 1/2 (or n = — 1), the coordinate transformations 

i^ l-^{l + bc)-b-fd\t-b-fip, 

f> = - l—^-lSlt + jip, 

z = z, (a -1/2), (55) 

bring the metric ([49| ) to the form 



2 



dsl = a^ e^" ('■+'') (dt^ - dr^) - dz^ -^{d^ + Sdtf, {a = 1/2). (56) 



a 



As shown above, this metric represents a flat spacetime but in a rotating coordinate system. 

It should be noted that the above coordinate transformations of Eqs.(p2h and (pq) are admissible only locally, if 
both sets of the coordinates are considered as representing cylindrical coordinates. Otherwise, it will give rise to new 
topology |2^J32]. However, in this paper we shall take the point of view that only the coordinates jx^^} — {t, r, z, ip\ 
represent the usual cylindrical coordinates. Then, the topological identifications discussed in |2l| , |32| are not applicable 
to the present case, so that the above coordinate transformations are admissible even globally. 



Assuming that a stationary shell located on the hypersurface r = tq > generates the spacetime described by the 
metric Eq.(|5J^) or Eq.(56), and that the spacetime inside the shell is vacuum and described by the metric Eg. ([48[) , 
where ro is a constant, we find that the first junction conditions that the metric are continuous across the surface, 
i.e., 



lead to 



for a 7^ 1/2, and 



9nu{r^ro) = g+i 



^A-^^ro), P ^ A-^<'^-'\ro), j ^ V^roA'^-\ro), S ^ n, {a y^ 1/2), 



Mro+d)^^^ 7 = V^^o, S^n, (a =1/2), 



for a — 1/2. Substituting Eqs.(58) and (^9|), respectively, into Eq. (p3|) and Eq. (jSq) , we find that 

dsl = B^''{r){de - dr^) - B^^^^'''^\r)dz^ - rlB^^^-^''\r)id^ + fldtf, (a ^ 1/2), 
for a ^ 1/2, where 



B{r) 



ro+d 



\l(2a-lf 



(^ ^ 1/2), 



and 



dsj^ = e 



2a(,-ro)(^^2 _ ^^2^ _ ^^2 _ ^2 (^^ ^ ^^^^2 ^ ^^ ^ -^y2), 



(57) 

(58) 
(59) 
(60) 

(61) 
(62) 



for a = 1/2. 

Now we are at the position of calculating the surface EMT of the shell. In the following, let us consider the two 
cases, a ^ 1/2 and a = 1/2, separately. 

Case a) o- 7^ 1/2: In this case, from Eqs.(|i|), (||), (|l]) and d^), we find that 



[X.r]- = 



2cr 



(2a-l)(ro + d)' 
(2a -1)2 (ro + d)' 



[W,r] = — ., ('7 7^1/2). 

ro+d 



(63) 



Substituting the above expressions into Eq.([2q), wc find 

d 



V^ 



ro {ro + d) ' 
9-0, 

4cr(l - a)ro - (2cr - ifd 
(2a-l)2ro(ro + d) ' 
4a2 



Pz 



Pip 



{2n-lY{ro + dy 



(a 7^ 1/2). 



(64) 



The combination of Eqs.(p4) and (B3^ shows that in the present case we have D > Q. Thus, now the fiuid is type /. 
For the shell to be physically acceptable, it has to satisfy some energy conditions, weak, dominant, or strong [g^ 
It can be shown that if the condition 



0<cr<l, d>0, 
holds, the weak energy condition will be satisfied, and if the condition 



(65) 
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d > \ f r'^ " "' (66) 

— I 4(7'' rp 1 ^ /T <• 1 

I (2,7-1)2' 3 < Cr S i, 



holds, the dominant energy condition will be satisfied, while if the condition 



io- rp 



0<^<i 



d> { (2'--!)'' "-" - 4' _ (67) 

holds, the strong energy condition will be satisfied. Clearly, for < a < 1, all the three energy conditions can be 
satisfied by properly choosing the free parameter d. It is interesting to note that in the case fi = it was shown that 
the solutions have physics only when < ct < 1/2 [^. Recently we extended it to the range a £ [0, 1] |^], which is 
exactly the same as that obtained above for the case with rotation, H. ^ 0. 

Case b) (T = 1/2: In this case, Combining Eqs.M, (Esl), (|6|) and (pO), we find that 



[V',r.]"-0, [x,r]"=«, M^-O, [T4^,,]---l, (^7-1/2). (68) 

Then, the surface EMT of the shell is given by Eq.(|^ with 

1 

^7= — , 
7-0 

9 = 0, 



Vz 



aro- 1 



p^ = a, (fT = l/2). (69) 

It can be shown that when a > the weak and strong energy conditions will be satisfied, while when < a < l/ro 
the dominant energy condition will be satisfied. Similar to the last case, the fluid now is also type J. 

V. CONCLUSIONS 

In this paper, the main properties of cylindrical spacetimes with rotation were studied by using the NP formulas. The 
physical interpretations of each component of the Riemann tensor was given by considering null geodesic deviations. 
It would be very interesting to see their experimental implications. 

In Sec. /// the spacetimes were extended to include rotating cylindrical shells, and the general expression of the 
surface energy-momentum tensor of the shells were given in terms of the discontinuation of the first derivatives of the 
metric coefficients. This would be very useful in studying gravitational collapse of a rotating cylindrical shell. As a 
matter of fact, this was one of our main motivations for such a study. As an application of the formulas developed, we 
considered a stationary shell that generates the Lewis solutions, which represent the most general vacuum solutions 
of the cylindrical spacetimes with rotation H, by assuming that inside the shell the spacetime is flat. It was shown 
that by properly choosing one of the free parameters, the shell can satisfy all the three energy conditions, provided 
< CT < 1, where a is related to the mass per unit length of the cylindrical shell (qP^. This range is exactly the same 
as that obtained in the static case [||. 

APPENDIX A: THE CHRISTOFFEL SYMBOLS AND THE EINSTEIN TENSOR 

Corresponding to the metric (Q), the non- vanishing Christoffel symbols, deflned by 

T^'a = -jS^" {.9tA,i. + gva,\ - gv\,a] , (A.l) 

are given by 



00 






^^'^"'' (^'* - -^) ' r?i = x,*-V',t, r?3 = -iw/^e-^x^,,. 
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^ 22 



o2(2^-x) 



V^,t, rl, = -w'e-'^U, 



w 



tL = Y.r - ^,r + uj^w\-^^ ( V;,, - -!^ - ^ ) , r^i = X,* - V',*, 






p3 

^03 



p3 
^ 13 



-{l~uj'W'e-'^) 



w 



^,i 



W.t 
W 



r?i = -w(x,t- V",*), 






^, 



)W^e'^^oj_r, r? 



-we 



2(2iA-x) 



v-,* 



r^g = cjw^2g-2x 






while the non-vanishing components of the Einstein tensor are given by 
Goo = 



^ - ^^X,W, + ^P,rW,r) + ^l + i^l - \oj'w'e-'^ul 



Goi — — 



Cj2y[/2g-2x 
W,tr 



X,tt - X,'. 



^-^^ 1.2 j,2 r,X,r'-^,T^ . o^.rW.r 



^^-fV^,-^,-2- 

UJ ■ ' UJ 



ojW 



1 /cj^y 

4 V w y 



1 



W ^^X,tW^r+X.rW^t) + 2i;,ti}^r 



1 



-iO 



w\-^^ 



UJ^tr — '^X,t^.r + +i^,r 



w^ 
w 



^03 

Gil 

Gi3 



Ti/2 -2x J , ^^rr , ,2 ,2 X,rW,r , 3uJ rWr 



2uM 



;W^2g-2x^2 









G33 = -VF^e-^x ^ ^_^^ _ ^^^^ + ^2 _ ^2^ ^ _M/2e-2x^2^ 



where x,r = dx/dr and x_4 = dx/dt, etc. 



(A.2) 



(A.3) 



APPENDIX B: THE SPIN COEFFICIENTS AND THE RICCI AND WEYL SCALARS 



Choosing the nuU tetrad as given by Eq.(0), we find that the spin coefficients are given by 



P = If^-^um^m" = - 
K = L-^m^T ^ 



c^-x 



2V2AW 



{W^t + W,r), 



a = l^-i^m'^ni'^ 



ei'-x 
2V2A 



2{'ljj,t+1p,r) 



Wt + W.r 



w 
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r = L-um'^n'' 



-i ^ LO . 



2V2 



a = - [l^j.-un^m" - m^^-i^m^m") = i 



We"^-^^ 



rt \''fi;iy"' "^ "t'jj-^i/iii' lit, J — o . nz ^•,'^1 



e = 2 ('^;'^"'''' ~ ™K-™'''") 



2V2A 



X,t + X,r - (-0,* + V'.r) 



^,i + Ar 



A 



t^^^ 



Ae'f'-^ 



/i = —nfj^-i^m^m 



V = —n^j^.i,m^n^ = 



TT — —n^.j^m^l^ — i 



2V2W 
I, 

Ae'f'-^ 

2V2 

2V2 



{W.t ~ W^r), 



2(V',t - 4>,r) 



W 



1 _ We^-"^^ 



1 _ Ae^^^ 



X,t - X.r - (V'.t - V'^r) + 



A 



Then, the corresponding Ricci and Weyl scalars are given, respectively, by 
1 



$ 



00 — T^'-'M' 



e2('/'-x) 
iA^W 



W^tt + 2W^tr + W,rr ^ 2(x,t + X,r){W,t + M^.r) + 2W {^ ,t + V',r)' 



1 



$01 = 2 V^""™' 



8A 



*02 = -5'^^m^m'' 



3w r 

W,rr + til^tr - 2uj^r{x,t + X,r) + 'T^i^-t + ^.r) 



^ ±e2(^-x) 



i^.tt - 0,rr + -^(t*M^,t - i^.rW^r) " ^(W^,« " W^rr) + ^W^e'^^LO^^. 



(B.l) 



$ii = ^^^.(^^n" + m^m'^) 



o2(V^-x) 



V',« - V',rr - i'l + ■4'1 + r-i^.tW^t " V-^rM^.r) + X,rr " Xm " -M^^g-^X^^^ 



$12 = -S^uU^m'^ 
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T4^ 



W,rr - W^tr + 2w,r(x,t " X,r) w"(^.* ~ ^''"^ 



$22 = 2'5m-^''«" 

^2g2(VJ-x) 



[VK,U - 2M^,t, + W^rr - 2{x.t - X,rKW.t - W.r) + 2VK(V',t " V,r)'] , 



A = — i? 
24 



12^ 



2(^-x) 



l/",** - V',rr - -0,4 + V',r + r^i'fp.tW^t - V'.rM^.r) " {X,tt " X,rr) 



W 4 



(B.2) 



and 



13 



where 



*o 



e2(V-x) 



*i 



2^2 

1 

" 2W 



{^p.tt + 2tP^tr + i^.rr + 2{tP,t + Ip^rf ~ 2{tP,t + 1p,r){x,t + X,r) 



[W,tt + 2W,tr + W,rr - 2{W,t + W,r)ix.t + X.r)] } , 



8A 



3w r 



UJ^rr + ^,tr - 2uj ^r{lp ,t + 1p ,r + X,t + X,r) + ~^i^,i + ^ -x) 



*2 = -\c^,uX5[l''n''l^n^ - l^n^'m^rh^] 



1 



e'('^~'^) {ij,tt - i'.rr - XM + X,rr + '^{tt " ^'r) 



+ ^ [T^,« - W.rr - Ki^,tW^t - i^,rW^r)] - W^e-'^^LO^ \ , 



A™<5 



*3 = -Cn^xsn^rn^fh 



i-AWe^'>'-^^ 



w. 



A^<5 



■ U,tr + 2uj^r{lp,t - Ip.r + X.t " X,r) 






iW.t^W,r) 



*4 = -Cf,^xsn''m''n^m 
1 



^2g2(^-x) 1^^^^ _ 2^,,,, + ^^,, + 2{tl;,t - ^,r)' - 2{i;,t - ^.r){x,t - X,r) 



1 

2W^ 



[T4^,„ - 2T^,tr + W,rr ^ 2{W,t - W^r)ix,t " X,r)] ^ , 



(B.3) 



and Ci_iucrS denotes the Weyl tensor. In terms of the Ricniann and Ricci tensors, it is given by 

Cfii^xs = Rfii^xs — -^ {gfixRi^s + 9i^sRfiX — guxRf^s — gtisRvx) 

-~-^{9iiS9ux - gf,xg^s)R- 

In terms of the Weyl and Ricci scalars, the Weyl and Ricci tensors are given, respectively, by 

Cf^^xs = -4 {(*2 + *2) {l[f,n^]l[xns] + m[^rh,,]m[xms]) 
- (^2 - *2) {l[p,n^]m[xms] + m[^fh^]l[xns]) 
+'^on[f,m^]n[xms] 

+*i {l[t^n^]'^[xms] + n[^rh^]l[xns] + nf^m^jTOfAm^j + TOf^m^jnfAm^]) 
-*2 {l[^m^]n[xms] + n[^m^]l[xms]) 

-*3 {l[t,nt,]l[xms] + l[f,m,,]l[xns] - ^[^m^jmfAm^] - m[^m^] /[Am^j) 
+-^4l[^m^]l[xms] 
+^oTHf,mi,]rnxms] 

+^1 {llt,n^]nixms] + nf^m^j/fAn^j + nf^m^jmfAm^] + mf^m^jnfAm^j) 
-*2 {l[f_cm^]n[xms] + n[^m^]l[xms]) 

-§3 (Z[^n^]Z[A?Ti5] + ll^Lmi,]l[xns] - ^[^m^jTOfATO^] - ?«[m'71^]^[a'ti5]) 
+'^4llt,m^]lixms]} , 

and 



(B.4) 



(B.5) 



(B.6) 
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(2$ii - 3A) {If^n^, + Un^) + (2$ii + 3A) (m^m^ + m^fhf,) 

-$12 (/pTO,, + Z,.?7i^) - $12 {l^^rriy + l„mf_,)} . (B.7) 



It should be noted that Eqs.(B.6) and (B.7) hold for the general case 
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FIGURE CAPTIONS 

Fig.l So and Sp are two infinitesimal 2-elenients spanned by 62 and 63 and orthogonal to the null geodesic C 
defined by l^^, passing So and Sp at the points O and P, respectively. A light beam meets So in the circle S. 




Fig.l 



17 



Fig. 2 (a) The image of the circle S on Sp is contracted because of the interaction of <I>oo- (b) The image of the 
circle S on Sp is deflected into an ellipse with its main major axis along 62 because of the interaction of ^o- 
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■«1 — ► e. 



(a) 




(b) 



Fig. 2 
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Fig. 3 A spherical ball consisting of photons cuts So in the circle S with the point O as its center. The image of 
the ball at the point P is turned into a spheroid with the main major axis along a line at 45'' with respect to ei in 
the plane spanned by ei and 63 because of the interaction of $1 and $01, while the rays are left undeflected in the 
e2-direction. 




Fig. 3 
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